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Abstract. In previous papers [1-2], we have presented hyperbolic governing equa- 
tions and jump conditions for barotropic fluid mixtures. Now we extend our results 
to the most general case of two-fluid conservative mixtures taking into account the 
entropies of components. We obtain governing equations for each component of the 
medium. This is not a system of conservation laws. Nevertheless, using Hamilton's 
principle we are able to obtain a complete set of Rankine-Hugoniot conditions. In 
particular, for the gas dynamics they coincide with classical jump conditions of con- 
servation of momentum and energy. For the two-fluid case, the jump relations do 
not involve the conservation of the total momentum and the total energy. 

Sommario. In precedenti lavori [1-2] sono state dedotte equazioni di governo iper- 
boliche e condizioni di salto per miscele fluide barotropiche. In questo lavoro tali 
risultati sono estesi al caso piu generale di miscele di due fiuidi conservative, tenendo 
conto delle entropie dei componenti. Si ottengono equazioni di governo per ciascun 
componente della miscela. Pur non essendo queste un sistema di leggi di conser- 
vazione, usando il principio di Hamilton si ottiene un insieme completo di condizioni 
di salto di Rankine-Hugoniot. In particolare per la gasdinamica, queste coincidono 
con le condizioni di salto classiche per la conservazione del momento e dell'energia. 
Nel caso dei due fiuidi, le condizioni di salto non coinvolgono la conservazione del 
momento e dell'energia totali. 
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1. Introduction 

Hamilton's principle is a well-known method to obtain the equations of motion in 
conservative fluid mechanics (see, for example, [3-6]). It is less known that this vari- 
ational method is suitable to obtain the Rankine-Hugoniot conditions through a 
surface of discontinuity [7-9]. The variations of Hamilton's action are constructed 
in terms of virtual motions of continua. The virtual motions may be defined both 
in Lagrangian and Eulerian coordinates [3,10]. Such virtual motions yield the gov- 
erning equations in different but equivalent forms. However the shock conditions 
are not equivalent. For example, for the gas dynamics, Hamilton's principle in the 
Lagrangian coordinates yields the conservation of momentum and energy. In the 
Eulerian coordinates, we obtain only the conservation of energy and the conserva- 
tion of the tangential part of the velocity field [7,8]. Here we use variations in the 
Lagrangian coordinates. 
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We assume that Hamilton's action is defined with the help of a Lagrangian 
function which is the difference between the kinetic energy and a potential depending 
on the densities, the entropies and the relative velocity of the components of the 
mixture. This potential can be interpreted as a Legendre transformation of the 
internal energy [1-2,11]. 

In [1-2,11], we have considered only the case of molecular mixtures. The hetero- 
geneous fluid, when each component occupies only a part of the mixture volume, 
can be described by the same system if one of the component is incompressible. 
Indeed, if the phase "1" is incompressible, the average density pi is related with 
the volume concentration of component (pi by pi = pwipi, where pio = const is 
the physical density of the phase " 1" . Hence, the knowledge of the average densities 
gives the volume concentrations and the physical densities. In the present paper, we 
do not distinguish these two cases. We shall call both cases "two- fluid mixtures". 

It is well known (see for example Stewart and Wendroff [12]) that the governing 
equations of two-fluid mixtures are not generally in a divergence form. In this case 
we may not obtain the shock conditions for the system. Moreover, the system is 
often non-hyperbolic, which means the ill-posedness of the Cauchy problem. The 
hyperbolic two-fluid models were constructed by many authors (see for example 
[13]). The problem to obtain the Rankine-Hugoniot conditions was an open question. 
This is the aim of our paper. By using Hamilton's principle in non-isentropic case we 
obtain the governing equations for each component and a complete set of Rankine- 
Hugoniot conditions generalizing those obtained in [1-2] for barotropic motions. To 
present the basic ideas, we consider first in section 2 the one- velocity case and extend 
this approach in sections 3 and 4 to the two-fluid mixtures. 

Let us use asterisk "*" to denote conjugate (or transpose) mappings or covec- 
tors (line vectors). For any vectors a, b we shall use the notation a*b for their 
scalar product (the line vector is multiplied by the column vector) and ab* for their 
tensor product (the column vector is multiplied by the line vector). The product of 
a mapping A by a vector a is denoted by A a. Notation b* A means covector 
c* defined by the rule c* = (A* b )*. The divergence of a linear transformation A 
is the covector divA such that, for any constant vector a, 

div(A) a = div (A a ). 

dz 

Let A be any linear mapping defined on fi () and B = —— be the Jacobian 
matrix associated with the change of variables z = M(Z), z belongs to SI. Then, 

div A = detB div ( — — a) , (1.1) 
V det B J 

where divo ( div ) means the divergence operator in SIq {SI). Equation (1.1) plays 
an important role. 

The identical transformation is denoted by 7, and the gradient line (column) 
operator by V (V*). For divergence and gradient operators in time-space we use 
respectively symbols Div and Grad. 

The elements of the matrix A are denoted be Oj where i means lines and j 
columns. The elements of the inverse matrix A^ 1 are denoted by a*. If f(A) is a 

df 

scalar function of A, the matrix — — is defined by the formula 

OA J 

fdiy = df_ 

\dAJi da)' 

The repeated latin indices imply summation. Index a = 1, 2 refers to the param- 
eters of components: densities p a , velocities u Q , etc. 
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2. One velocity fluid 



The consequences of this section are well known. We obtain the classical governing 
equations and the Rankine-Hugoniot conditions for the gas dynamics. Nevertheless, 
the presented method is universal and is extended for two-fluid mixtures in the 
following sections. 

Let z = (^^j De Eulerian coordinates of a particle and T>(t) a volume of the 

physical space occupied by a fluid at time t. When t belongs to the finite interval 
[io,fi], D(t) generates a four-dimensional domain £1 in the time-space. A particle 
is labelled by its position X in a reference space T>o- For example, if T>(t) consists 
always of the same particles V — T>(to) and we can define the motion of the 
continuum as a diffeomorphism from T>(to) into T>(t): 

x = ¥Jt(X). (2.1) 

The motion (2.1) of the fluid is generalized in the following parametric form 



t = g(X,X) 
x = 0(A,X) 



M(Z) 



(2.2) 



where Z = ^ ^ ^ belongs to a reference space denoted Qo and M is a diffeomor- 
phism from a reference space Qo into the time-space ft occupied by the medium. 
Equations (2.2) lead to the following expressions for the differentials dt and dx, 



(£)-»(£) 



where 



r r 2i 

dZ dj>_ 
. dX 

In an explicit form, we obtain from (2.3)-(2.4) 

dt = 



dg_ 
<9X 

d<j> 

<9X J 



(2.3) 



(2.4) 



(2.5) 



Eliminating d\ from the first equation of (2.5) and substituting into the second, 
we obtain 

dx = u dt + FdX, 
where the velocity u and the deformation gradient F axe defined by 

d(f> ,dgx-i d<j> d<j> dg,dg^-i 



Let 



dX V <9A 



t = G(A,X,e) 



<9X 



or 



dX dX^dX' 



z = M £ (Z), 



(2.7) 



[ x = $(A,X,e) 

where e is a scalar defined in the vicinity of zero, be a one-parameter family of 



virtual motions of the medium such that 

M (Z) = 



M(Z). 
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We define the virtual displacement £ = (r, £) associated with the virtual motion 
(2.7): 

dG 9$ 
r= ^-(A,X,0),€ = ^(A,X,0) 

or, C =^(Z)|^o 

From the mathematical point of view, spaces fio and f2 play a symmetric role. 
From the physical point of view they are not symmetric: the tensorial quantities 
(thermodynamic or mechanical) are defined either on Oo or on Q. Their image in 
the dual space depends on the motion of the medium. For example, the potential of 
body forces If is defined on Q and the entropy s is defined on Qq. 
Let us consider any tensorial quantity represented in the form 

(t,x) e n — /(t,x). 

The tensorial quantity associated with the varied motion is 

/ (A,X, £ ) = / (G(A,X,e) , *(A,X,e)) = /(M e (Z)). 
We define the variation of / by 

5f= ^ (A,X,0). 

For a tensorial quantity represented in the form 

(A,X) G Q Q —> ft(A,X), 
the tensorial quantity associated with the varied motion is unchanged and 

Sh = 

Let the Lagrangian of the medium be defined in the form 

L = L(z, ^ , z) =L(z,B,Z) . 

This expression contains the gas dynamics model where the Lagrangian is [3] 

L= i p(l + |u| 2 ) -s(p,s)-pU(z) 

V ' . (2.9) 

= -pV*v- e ( P , s )- P n(z) 

Here V = i s tue time-space velocity, II(z) is the external force potential, p 

is the density defined by 

p dot F — p (X), 
and s is the entropy per unit mass defined by 

s = s (Z). 

It is not necessary to assume that so is a function of X only. This property will 
be a consequence of the variational principle (see formula (D.l) in Appendix D). 
The Hamilton action is: 

L(z,B,Z)dQ. (2.10) 
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For the gas dynamics we obtain from (2.4), (2.6): 

v-(i)-£ 

where t = (1,0,0,0) and fj, = t Bt = b\= |f. 

Consequently, 

2 ( 1 + l O = 2 ^ 

Moreover, 

P = ^ P0(X). (2.13) 
In the Lagrangian coordinates Hamilton's action (2.10) is 

L(z,S,Z) detB dfi , 

and the varied action is 

•M-j^M.M, ^,z)de«(^),BV 

Let T(n) be the tangent bundle of £2. 
The Hamilton principle is: 

For any continuous virtual displacement £ belonging to T(SY) such that £ = 
on T(dfi), 

Sa= — \ e = o = 0. 

Consequently, 

<Sa = / jdetB |^C + detB tr( + L <5(det B)}dfi . 

The Euler-Jacobi identity 

SdetB = tr( dA ^ B SB^J =tr (B^ 1 det B8B) = det Btr(B^ 1 SB) 
and the relation issued from definitions (2.4), (2.8) 



yield 



where 

S* = § and T = LI + BH. 
The Gauss-Ostrogradskii formula involves 

8a = / (S* - Div T)C d fi + / N*TC dej, 

Jen 

where N* is the external normal to dfl and dtu is the local measure of <9fl 
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If the motion is continuous on Q and £ = on dQ , wc got 
5a = / (S* - Div T)C d CI 



Consequently, the governing equations are 

S* - divT = 0. (2.14) 

In Appendix A we verify that (2.14) corresponds to classical momentum and energy 
equations. If there exists a surface E of discontinuity of B separating Q into two 
parts fli and C7 2 we get 



8a 



I (S* - DivT) Cdfii + / (S* - DivT) (dft + / N*[T]Cdu;, 



where [T] = Ti - T 2 denotes the jump of T S. 

Consequently the fundamental lemma of calculus of variations involves the Rankine- 
Hugoniot conditions 

N*[T] = 0. (2.15) 

Because N* is collinear to [— D n , n*], where D n is the normal velocity of E 
and n is the normal unit space vector, relations (2.15) are the classical Rankine- 
Hugoniot conditions representing the conservation of momentum and energy through 
the shock (see Appendix A and [7]). 



3. Two-fluid models: General calculations 

We shall study now two-fluid motions, the method being extended to any number 
of components. We generalize the representation of the motion (2.2) considering the 
motion of a two-fluid mixture as two diffeomorphisms [14] 

z = M a (Z a ), 
belongs to a reference space Q a associated with the a— th 



where Z a = 



component. The Jacobian matrix is defined by the formula 

Ba - az7 (Zq) - 

The velocity u a and the deformation gradient F a are defined similarly to (2.6). 
Two one-parameter families of virtual motions are associated with the two diffeo- 
morphisms 

(z = M Ml (Zi) r , n 

\z = M 2 (Z a ) 

such that Mi, (Zi) = Mi (Zi), and 

r z = Mi (zo r321 

\z = M 2 , £2 (Z 2 ) 
such that M 2 , (Z 2 ) = M 2 (Z 2 ). 

The two families extend the concept of virtual motion defined in section 2. Consider 
the family (3.1), but all the consequences are the same for the family (3.2). We 
define the virtual displacement of the first component by the relation 

9Mi, ei 

Ci - -a^-(Zi) k=o- 
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Virtual motion (3.1) generates a displacement of component "2" 

Mi, £l (Zi) = M 2 (Z 2 ), 

which defines Z2 as a function of Zi and si. Taking the derivative with respect 
to si and denoting 

r 7 dZ 2 

SlZ2 = ai7 k=0 ' 

we obtain 

5iZ 2 = B^d- (3.3) 
Let us consider any tensorial quantity / in Eulerian coordinates: 

z G Q —> /(z). 

The tensorial quantity associated with the varied motion is then 



/(Zi.ei) = /(M Ml (Z0) 



and consequently Sif = -^-(Zi,0) is the variation of /. 

OEl 

Let us consider the Lagrangian of the medium in the representation 

L = L(z,B 1 ,B 2 ,Z 1 ,Z 2 ). (3.4) 
For example, the Lagrangian of a two-fluid mixture is [1-2, 4-6, 11]: 

L = \ pi (l+|ui| 2 ) + i P 2 (l + |u 2 | 2 ) - W(pi,pa,«i,«2,U2-ui) - pn(z), (3.5) 

where p a is the density of the a — th component defined by 

p a det F a = p ,a(X a ), 
s a is the entropy per unit mass of the a — th component defined by the relation 

Sa — SO,a(Z a ), 

and n(z) is the potential of external forces. 

One can easily obtain the formulae analogous to (2.11)-(2.13) for p a and u a in 
terms of B a and Z a . Hence, the Lagrangian (3.5) may be rewritten in the form 
(3.4). The variation associated with the application (3.1) yields 

S ia = / Si (Ldet-Bi) dQi. 

Calculations presented in Appendix B give the following result: 

4 '° = /„{' r ( T ^- T '^) + s: <'}* (3 ' 6) 



where 



Ti = B 2 ^-B 2 , (3.7) 

U-D2 

_ dL , dL 



Meccanica_extended.tex; 11/02/2008; 1:14; no v.; p. 7 



8 Henri Gouin and Sergey Gavrilyuk 

The Gauss-Ostrogradskii formula and relation (3.3) involve 

5ia = \ St-Div T + (Div T^B^ 1 I Ci dfl+ / N*(T-TiB 2 _1 ) d dE. (3.8) 
./n J Jan 

Using the arguments described in section 2, we obtain from (3.8) both governing 
equations and Rankine-Hugoniot conditions for component "1" 



SI — Div T + (Div T^B^ 1 = 0. 
N*[r-TiB 2 _1 ] = 0. 
Since T — TiB^ 1 = L I + B\ -J^-, equation (3.10) is equivalent to 



N* 



LI + Bj 



dL 



= o. 



(3.9) 
(3.10) 

(3.11) 



Equations for component " 2" are obtained by permutation indexes " 1" and " 2" : 

S* 2 — Div T + (Div T 2 )Br 1 = (3.9') 

and 



N*[T-T 2 Sr 1 ] = 
Formula (3.10') can be rewritten in an equivalent form 

dL 



N* 



L I + B 2 



8B 2 



= 0. 



(3.10') 



(3.11') 



Let us remark that the equations (3.9) and (3.9') are not in a divergence form. Let 
us denote 

Q » dL 



and 



Identities 



So Q — det B a — , To a = — det B a B a 

OZj a Cilia 



(3.12) 
(3.13) 



T 02 = -(detB 2 )B 2 x Ti , T i = -(detBi)^ V2, 

1 1 

S i = S * + dclB^ S( ' 2jB2_1 ' = S * + dclBT S ° 1jBi " 
and (1.1) involve that (3.9) and (3.9') are equivalent to 

S*-DivT+ _J_ (SSi-DiviToi)Bf 1 = 0, 
det ±>i 



and 



S* - Div T+ — I— (S,* 2 - Div 2 T m )B^ = 0, 
aet _£>2 



(3.14) 



(3.14') 



where Div a Tb a means the divergence of To a with respect to the a— th Lagrangian 
coordinates. The following identity (C.l) proved in Appendix C, 

S* — Div T + — (SSi - Divi ToOBf 1 + 775- (SS2 - Div 2 T 02 )S 2 - 1 = 
det Si deti?2 



and (3.14), (3.14') yield 



S* - Div T = 0. 



(3.15) 
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We see in the next section that for fluid mixtures, equation (3.15) represents the 
conservation laws of total momentum and total energy. Hence, (3.14) and (3.14') 
are equivalent to: 

Soi -DiviToi = 0, (3.16) 

Sq 2 -Div 2 T 2 = 0. (3.16') 

Equations (3.16) and (3.16') represent equations of motion for each component of 
a two-fluid medium in Lagrangian coordinates. They are less useful than equations 
in the Eulerian coordinates. However, they are in a divergence form and involve the 
conservation of the total momentum and the total energy. 



4. Application to two-fluid mixtures 



For a two-fluid mixture, the Lagrangian is (see (3.5) 



L = ipilVrl 2 + ip 2 |V 2 | 2 -M/( Pl ,p 2 , Sl , S2 ,w)- P n(z), (4.1) 



where 



V = (M, *lt, (°) = V 2 - Vl 



B a 

Ha 

pa pOa (X a ) 



(4.2) 



(4.3) 
(4.4) 

Definitions (4.1)-(4.4) involve the governing equations (3.9), (3.9') in the following 
form (see Appendix D): 



det B a 

Set = S() a (Z a ). 



9Ka 

dt 



+ rotK a x u a = V*(R a - K* u a ) + 6 a V*s a , 



^ + div (p a u a ) = 0, 



(4.5) 



d{p a s a ) 
dt 



+ div (p a s a u a ) = 0, 



where 



p a K a = - — =p a u a -(-l) — — , with w = u 2 -ui, 
au a ow 

1 2 dW 

R a = ^ — = - u„ - n, 

Op a 2 Op a 

a dL dW 
p a a a = — — — — — — . 

os a os a 

We notice here that the governing equations were obtained earlier by a different 
method in [11]. Equations (4.5) yield the conservation laws for the total momentum 
and the total energy associated with (3.15) 

^(piKi +/92 K|)+ 

/ dW dW \ 

+div (pniiK? + p 2 u 2 K; + (pi — + p 2 — - W) I J + pVn = 0, 
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a=l 
2 



+div P«u a (K* u Q - - p — = 



(4.6) 



In the general case, they are the only conservation laws admitted by the system. 
Hence, the system (4.5) is not in a divergence form. The Rankine -Hugoniot condi- 
tions (3.10)-(3.10') (or (3.11)-(3.11')) for this system are obtained in Appendix A(b) 
(see formulae (A. 9)); they are 



- D n (L - p a K*u a ) + p a n* u a (R a - K*U Q ) = 0, 



- D n p a K a + n* ((L - p a R a ) I + k;) 
In addition, the mass conservations laws are expressed in the form 



0. 



p Q (n*u Q - D„) 



= 0. 



(4.7 1 ) 



(4.7 2 ) 



(4.7 3 ) 



Let us consider shock waves when n*u Q — D n ^ 0. Taking into account (4.7 2 ) 
and (4.7 3 ) we obtain 

(4.8 1 ) 



0, 



Ka - (K Q n)n 

which means that [K a ] is normal to the shock. By using (4.7 1 ) and the identity 

1 



p a u a (R a — K*u Q 



L — p a K* u a 
we obtain from relation (4.7 2 ) 

K*u a - R a - D n (K*n) 
Consequently, (4.7 1 ) yields to 

L - p a Ra + pa(n*U a - D n )K a U 



0. 



(4.8 2 ) 



(4.8 3 ) 



Equations (4.7 3 ), (4.8 1 ) - (4.8 3 ) form a complete set of eight scalar Rankine-Hugoniot 
conditions, representing the conservation of mass, momentum and energy of a- 
th component. For the gas dynamics, these conditions correspond to the classical 
shock conditions. We have to emphasize here that for the two-fluid model equations 
(4.8 1 ) — (4.8 3 ) do not involve the conservation of the total momentum and the total 
energy through the shock. We notice that the jump conditions (4.8 1 ) — (4.8 2 ) were 
obtained earlier for barotropic fluids using the variations in the Eulerian coordinates 

For contact discontinuities, when n*u Q — D„ = 0, we get from (4.7 ) — (4.7 ) 



L — p a Ra 



0. 



(4.8 4 ) 
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For the gas dynamics equation (4.8 4 ) corresponds to the continuity of the pres- 
sure. All the jump conditions are obtained from the Hamilton principle without any 
ambiguity. 



5. Conclusion 

Using the Hamilton principle we have obtained in the general conservative case 
the governing equations for two-fluid mixtures (4.5). The equations for the total 
quantities ( total momentum and total energy) are in a divergence form (see (4.6)). 
The equations for the components are in divergence form only in the Lagrangian 
coordinates (see (3.16)-(3.16')). The Hamilton principle gives also a set of Rankine- 
Hugoniot conditions (4.7 1 ) — (4.7 2 ). Together with the equations of mass (4.7 3 ) they 
form a complete set of the jump relations. For the gas dynamics model they coincide 
with the classical Rankine-Hugoniot conditions of conservation of mass, momentum 
and energy. 



6. Appendix A 

In gas dynamics, the Lagrangian (2.9) is a function of V = [ ] , p, s and z. 



, u j 

Formulae (2.11)-(2.13) allow us to consider the Lagrangian as a function of z, B, Z: 

V = — £ , p = £*B£, 

M . (A.l) 

p det B = p po(X), s = so(Z) 

We need to calculate T = L I + B . From (A.l) we have: 

dB 

dp = t dB £ = tr(£t dB), (A.2) 

dV = — t- ^ dp = — £- ^tr (WdB), (A3) 
p p 2 p p 2 



£ = ^^(^j+MX)^^ = p { —-B (AA) 



Moreover, 

ds 

dB = °" 

The differential of L = L(V,p,s,z) is: 

dL „ r dL , dL , dL , 
dL = — dV + — dp+— ds + — dz. 
ov dp ds dz 

Taking into account relations(A.2)-(A.5), we obtain 

dL f dB £ B£ 



dL , dL , dL , 
+ dp- dp+ ds- ds+ dz- dZ = 



(A.5) 
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jidW 



dL 
dp 



A 4 

dL dp dL ds 



1 TvTT^ ' 77777z.J' /Z ~ f/z " 



dL 
to' 



Hence, 



Then 



dL 
dB 



£ j£r U "(^ m ) dLdp_ 
+ dpdB' 



'aw 



(A6) 



, dL 



B dB- 



/' 



Taking into account (A.l), we obtain 
„dL 



<9B dW \ dW J ^ dp 

and, (A. 7) yields then 



<9l 



vr. 



(A.8) 



i 



a) In the case L = -pV*V — s(p,s) — pll formula (A.8) can be rewritten in the 
form 

T = pi + P ww* - Qv*v + 4 + n)pvr, 

where p = pej, — e. Then, 

—E , pu* 
-(S + p)u , p I + puu* _ 

where E = ^pV*V + e + pll is the total volume energy. 

Moreover, since S* = — p^j— , we get the energy and momentum equations for gas 
dynamics motions 

'§+*(<^>»)-f = °' 

dpu* ,. , T 9n 

+ div (p 7 + puu* + p— = 0. 
k at ox 

Because N* is collinear to ( —D n , n* ), the classical Rankine-Hugoniot conditions 
(2.15) are 

D„E- (E+p)n*u = 0, 

D n p u* — n* (p I + puu 
b) For the two-fluid model the Lagrangian is 

L = ipiV*Vi + ip2V*V 2 - W(pi,p 2 ,si,s 2 ,w) 
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with 



(°) = ( ° )=v a -v 1 

\W J \U2 - Ul j 

We obtain similarly to (A. 8) 



Here 



dB a 

dL 
dV a 



dL_ _ ( dL_ 

+ a av a Uv Q 



dL dL 



a pen 



» aw * aw aw _ * 



av D 



9w <9V Q 



ax _ i v * v aw 



In a matrix form, we have 
dL 



LI + B a - 



dB a 



L — /9 Q K* u Q , p a K* 

"paU a K a U a -(- PaRa U a , (L — p a Ra) I + p a U a K* 

and the Rankine-Hugoniot conditions for the two-fluid model are 

L — paK^Ua , paK*a 



N* 



p a U a (R a - K* U a ) , (L - p a R a ) I + PaU a K a 



= 0. 



(A9) 



7. Appendix B 

The variation 8\a is: 



5 ia = J j^iL + L <5i(detBi)(detBi)" 1 |dfi = J ^L + tr^L^- 



Si 



Formulae (B.3) and (B.4) involve that 



>dn. 



Formula (3.3) involves 

i- t dL . „ dL . „ dL aL ^ , „. 

W = *r|^a + gg^ft J + _B 2 ^ + ^C, (B.l) 

By definition, 

SiB! = ^Bi. (5.2) 

Let us calculate <5i-B 2 . We get 

-* ( &) -* -* (&) (H) • <-> 

Since for any linear mapping A : 5(^4 ) = -A^SA A~ , we have 



azi \ _ azi /dz 2 \dZj dz 1 ds 1 z 2 dz 1 

az 2 ) az 2 1 v azi / az 2 az 2 azi az 2 ' [ ' ' 



ar-^ir a z a^iZa _ ad R „ a^z 2 D , . 
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Substituting (B.2) and (B.5) into (B.l) we get 

r T t ( dL d<, „ dL ( dC 9<5iZ 2 \\ dL dL 

-fr/fR 9L +R 9 M 9Cl R 9L B d& ^\+ 9L R-V + 97 V 
-' r jl Bl 9BT +B2 -fr- B2 dm B2 ~dz~ f + dZ 2 B * Cl+ fe Cl " 

Hence, 

where T , Ti and Sj are defined in (3.7). The relation (3.6) is proved. 



8. Appendix C 



Theorem: The following expression is an identity 
S* -Div T + ( S ' 1 -DivToi) Bf 1 + (S^ -Div T 02 ) Bj 1 = 0, (C.l) 

w/iere (see formulae (3.6), (3.12) -(3.13)) 



S* 



SS a ^detB^ 



9Z Q 



To a = — det B a 



dL 

dB a 



B a . 



Proof: Using (1.1), we get 



CK — 1 CK — 1 

We have to prove that this expression is identical to 

n ,r- f . n ,( I , + .,| + ^)-« 



Indeed, by definition 

&L 9L aZi <9L <9Z 2 



_ j r dL dL OZi #L yz 2 ( dL dB 1 \ / 
GradL=— + —— + —— +tr — - — - +ir 
c*z 9Zi dz 5Z 2 5z V oBi <9z / V 



dL dB 2 
dB 2 dz 



where 



tr 



dL dB a 



dB a dz 

Hence, we have to verify only the formula 



dL V dB]_ 
dzP ' 



dB a 



(C.2) 
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dz 

To prove the identity (C.2) we may consider only L — L(B), with B — -jr=. 

dZ 

dL , 
Let A — ——, then 
oB 

dL . ( dL dB \ db] 



\ dB dz m 



dz m "' V dB dz m I aj dz-" 



Hence, 

d{b) 4 b k a ) - _ d{b) a{) k t , db k _ 

d? p ~ d? 3 p + j k d7 p ~ 

d ( h ) «p) , h i j db k s dZ m - s _ 9(6} a p ) . . db k s r 



dz- + 3 k dz^ ~~d^ v = dz- + 3 1 °* dz^ p = 

4) . m db k m - _ d{b) a p ) m db k m ,- s _ d{b) a p ) m db k m 
—q— +a k Q^b p - + a k -Qjb s b p - + a k — 

This is a proof of (C.2) and consequently (C.l). 



9. Appendix D 



First of all, we obtain the governing equations for each component in the Lagrangian 
coordinates. These equations yield easily the governing equation in the Eulerian 
coordinates. In the Lagrangian coordinates equations (3.16) — (3.16 ) are 



S<j a - Div a T * a = 0, 



where 



and 



c* -i . 75 dL dp 0a 

dZ a dZ a 



detB a p a 8, 



ds Qa 
dZ a 



dL 



Toa = — detBa, — B a . 

dB a 



Similarly to (A. 6) we obtain 



dL P a ^lW^ 



dB a p a 
Consequently, 



It . dL . 
p a 9V„ 



It 

Pa 



Br, 



detB « TTEr Ba = P ()a l itrr Ba ~ p°<* (Jr v o)^* B <" + po a R a {M*B a - p a /). 

Ct-Da C V a O V a 



If \ a — t, p a — 1, B a 



1 
II,, F a 



Toa — POa 



and consequently 

, -k;f c " 

, R a I 



The governing equations of a-th component in the Lagrangian coordinates (t, X a ) 
are: 

-p a da dct B a = 0, 

^(p 0a K.* a Fa) - AlVa(pOaRa I) + Ra - POaOa = 0. 
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Therefore, 

^=o, <zu> 

^(K a F a ) - V a R a - e a V a soc = 0. (D.2) 
Taking into account the identity 

dF a _ du a 

dt ~ 9x a 

and substituting the partial derivative with respect to time in the Lagrangian coor- 
dinates by the material derivative in the Eulerian coordinates, we obtain 

|. a = o, with | = | + Ua -v, ( mo 
^ k; + k; ^ = vR a + o a vs a . (d.2') 

Let us note that (D.l'), (D.2') and the mass conservation laws are equivalent to 
(4.5). 
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